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Abstract: We experimentally investigate quadrature correlations between 
pump, signal, and idler fields in an above-threshold optical parametric os- 
cillator. We observe new quantum correlations among the pump and signal 
or idler beams, as well as among the pump and a combined quadrature of 
signal and idler beams. A further investigation of unforeseen classical noise 
observed in this system is presented, which hinders the observation of the 
recently predicted tripartite entanglement. In spite of this noise, current 
results approach the limit required to demonstrate three-color entanglement. 
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1. Introduction 

The optical parametric oscillator (OPO) is a very well-known source of non-classical light. 
Quantum correlations between the intensities of the downconverted beams, signal and idler, 
were first measured already in 1987 [1]. Because of the high degree of correlation attain- 
able between these two bright light fields they were called, in the above-threshold opera- 
tion, twin beams. In fact, squeezing in the intensity difference reached the impressive value 
of —9.7 dB [2]. As was first recognized by Reid and Drummond in 1988, the quantum corre- 
lations are not restricted to intensity but also extend to their phases; the OPO produces twin 
beams in an Einstein-Podolsky-Rosen (EPR) state [3]. EPR-type entanglement was obtained 
a few years later in the OPO, but only when operating below threshold [4]. Since then, this 
system has been used in many applications in continuous-variable quantum information [5-7]. 

The measurement of entanglement in the above-threshold OPO would take 17 years to be 
realized [8-10], in part because of technical difficulties to perform phase measurements on 
bright beams that usually have different frequencies. Another important reason was the presence 
of an unexpected excess noise generated in the intracavity pump field by the non-linear crystal, 
which was transferred to the sum of downconverted phases during the parametric process [11]. 
As a consequence, squeezing in this observable could only be found very close to the oscillation 
threshold [8]. This strange feature, unaccounted for by the theoretical model, went unnoticed 
until recently, when phase measurements were performed above threshold. We present here 
a further study of the fields' quadrature noises as functions of pump power to show that the 
standard theoretical model for the OPO can not account for the observed features. Even the ad 
hoc model presented in ref. [11] reveals itself to be insufficient. These noise features warrant 
investigation in themselves, which will have to go beyond the present work. 

Another "surprise" was to be encountered in this thoroughly studied system. Although 
quadrature squeezing of the reflected pump beam was experimentally observed in 1 997 [12], lit- 
tle attention was paid to quantum properties of the pump field. In many theoretical approaches, 
the pump was treated as a classical quantity. A very recent investigation of the quantum proper- 
ties of the full three-mode system, by Villar et al. revealed that all three fields should be entan- 
gled [13]. This seems to be the simplest system predicted to directly generate genuine tripartite 
continuous-variable entanglement. It also appears as a distinctive feature of above-threshold 
operation, in contrast to the widely-used OPO below threshold. The first observation of a triple 
quantum correlation was recently reported by Cassemiro et al., who observed squeezing in 
the combination of the sum of twins' phases with the reflected pump amplitude, for a detuned 
OPO [14]. This result was a good indication that phase-phase correlations should be present 
for a resonant OPO, as theoretically predicted [13]. 



Here we present the first direct observation of such triple quantum phase correlations. For 
higher pump powers, we show as well that the amplitude of the reflected pump beam is 
quantum-mechanically correlated to each of the twin's amplitudes. As for the tripartite en- 
tanglement, owing to the spurious excess phase noise introduced by the non-linear crystal in 
the intracavity pump field [11], the quantum correlations observed do not suffice. 

The quest for bright multicolor entangled beams will be rewarded by providing an interesting 
tool for quantum information networks, since they can concomitantly interact with physical 
systems otherwise incompatible. We describe below our latest experimental advancements in 
this direction. 



2. Tripartite entanglement 

The existence of direct tripartite pump-signal-idler entanglement can be physically understood 
as a consequence of energy conservation. The parametric process imposes for the frequencies 
of these fields the relation Vo=ViH-V2 (indices 0, 1, and 2 refer to pump, signal, and idler fields, 
respectively). Phase fluctuations, regarded as small frequency fluctuations, are therefore con- 
nected: 5(po=5(pi+8(p2. This leads to strong quantum correlations among the phase fluctuations 
of the three fields. On the other hand, the creation of a pair of signal and idler photons inside 
the OPO cavity only occurs upon annihilation of a pump photon. In triply resonant cavities, as 
we consider here, pump depletion above threshold is a well known effect. Correlations between 
the amplitudes of the pump and twin beams are thus created. 

In order to demonstrate genuine multipartite entanglement, a more rigorous approach makes 
use of inequalities which must be fulfilled by all separable states. A violation of a subset of the 
inequalities suffices. We denote the electromagnetic field's amplitude and phase quadratures as 
Pj = exp{—i(pj)aj +exp{i(pj)d^j and qj ~ -~i[exp{—i(pj)aj — exp{i(pj)dj], respectively, where 
j G {0, 1,2}. The fields' annihilation and creation operators are denoted aj and a'j, satisfying 

= 5j f, f E {0, 1,2}, and (pj are the fields' mean phases. Tripartite entanglement among 
three fields can be witnessed by violation of the following inequalities [13, 15]: 

where aj are numbers that minimize Vj. These are generalizations for tripartite systems of a 
criterion for bipartite systems derived independently by Duan et al. and Simon [16, 17]. In fact, 
if one makes aj = 0, the bipartite criterion is recovered. Violation of at least two inequalities en- 
tails genuine tripartite entanglement. In our case, one has V\ ~ V2 according to the model, since 
signal and idler beams play the same role in the OPO equations, differing only by polarization 
for type-II conversion, and by optical frequency. In the experiment, these two inequalities dif- 
fer slightly owing to imperfections, such as unbalanced losses for signal and idler intracavity 
modes. 



The first inequality contains well-known correlations. Its first term, A /3_ [where /3_ — {p 



P2)/V2] corresponds to the amplitude correlation between the twin beams. The second one, 
A^q'_^_ [where q'^ = {qi +q2)/\/2 — Oloqo]^ can be written as the difference of A^q+ [where 
q+ = (^1 + q2} / V2], which involves the twins' phase anti-correlations, and a correction term 



j3o given by their phase correlation with pump, 



A q_^_ -A q+- Po , Po lA^qo ' 

where Q^^y = {Sqjdqf), f ^ j. Notice that j3o > 0. If A^q'^ < and j3o ^ 0, there is a quantum 
correlation between all three fields. The inequality corresponding to Vq is violated even in the 
absence of the correction term, since signal and idler are entangled [8]. The correlation to the 
pump, signaled by )3o, strengthens this entanglement and increases the region of parameters 
where violation occurs [14]. 

Here, we present an extended experimental investigation of this system, by measuring in ad- 
dition the quadrature noises appearing in the second and third inequalities. To simplify notation, 
we named the terms appearing in Vi and V2 as 



Following Eq. (IH, we write 



where ^o; = i^j — Qo) /V2 and / 7^ j. In this case, qoj is the composite noise of the pump and 
one of the twin beams, and j3y/ is a correction term directly dependent on their correlations with 
the remaining beam. 

3. Experimental set-up and system characteristics 

The optical parametric oscillator comprises a non-linear j'^' crystal disposed inside a linear 
Fabry-Perot cavity. The coupling mirror for the pump beam is a partial reflector at 532 nm 
(reflectivity R =69.4%) and a high reflector (R >99.8%) at 1064 nm; the coupling mirror for 
the twin beams has reflectivities which are equal to 96.0% at 1064 nm and greater than 99.8% 
at 532 nm. The crystal is a 12 mm long type-II High Gray Tracking Resistant Potassium Titanyl 
Phosphate (KTP) from Raicol Crystals Ltd. The crystal temperature is kept constant by a peltier 
element near 25 °C, with a stability of the order of tens of mK. Its losses are approximately 
3% at 532nm and less than 1% at 1064nm. Typical threshold power is about 75 mW, and the 
OPO cavity bandwidth is 5v„p„=45(2) MHz. A frequency doubled diode-pumped Nd:YAG 
source (Innolight Diabolo) at 532 nm is used to pump the OPO. This laser is 'filtered' by a 
mode-cleaning cavity [18] (bandwidth of 2.4 MHz) prior to injection in the OPO, in order to 
remove all classical noise which could hinder the entanglement for analysis frequencies above 
15 MHz [11,19]. 

Phase fluctuations are measured through a self-homodyne technique using optical cavi- 
ties [19, 20]. An optical cavity adds a frequency-dependent phase to the field it reflects, as a 
function of the detuning between the incident beam carrier frequency and the cavity resonance 
frequency. Scanning the cavity resonance (e.g. with a piezo-electric device to change its length) 
results in different dephasings between carrier and sidebands for each cavity detuning. In the 
Fresnel plane this can be qualitatively visualized as a rotation of the noise ellipse with respect 
to the mean field. As a result, the intensity fluctuations of the beam reflected by the cavity 
can correspond to fluctuations of any desired quadrature of the incident beam. It can be shown 



that the incident phase quadrature is totally converted into the reflected amplitude quadrature 
when the analysis frequency Q. is higher than \/2 times the cavity bandwidth 5v [19,20]. The 
cavity-carrier detunings A (relative to the cavity bandwidth) for which it occurs are A w ±0.5 
andAw±n/5v. 



.Analysis Cavity' 1 



OPO 




rP'^ Rotator i J7 T" 

^ V I TA KTP zJ 




to demodulating 
chain... 




.Analysis Cavity 



Fig. 1. Sketch of the experimental setup. PBS: polarizing beam splitter. 

A sketch of the experimental set-up is presented in Fig. [T] There are three analysis cavities 
in the apparatus, with the following bandwidths: 5vo=l 1.5(6) MHz, 5vi=14.5(l) MHz and 
5v2= 13.6(2) MHz, where the indices correspond to the beam analyzed in each one. In order to 
measure the pump beam reflected by the OPO we use a Faraday Rotator and a polarizing beam 
splitter (PBS). The infrared beams are detected with high quantum efficiency (93%) photodi- 
odes (Epitaxx ETX300), resulting in an overall detection efficiency of rj = 87(3)%. Another 
photodetector is used for the pump beam (Hamamatsu S5973-02, quantum efficiency of 94%), 
with an overall detection efficiency of 770 = 74(3)%. 

The photocurrent high frequency component (HF) is sent to a demodulation chain, where it 
is mixed with a sinusoidal reference at the desired analysis frequency (0 = 21 MHz), within a 
600 kHz bandwidth. Analysis frequency is chosen to have the smallest value (since the quantum 
correlations increase for decreasing frequency) such that the analysis cavity can still completely 
rotate the beams' noise ellipses. The resulting low frequency beat note is acquired at a rate of 
600 kHz by an analog-to-digital (A/D) card connected to a computer Variances of each indi- 
vidual noise component are then calculated taking a group of 1000 acquisition points. Finally, 
the noise power is normalized to the standard quantum limit (SQL), previously calibrated [14]. 



4. Measurement results 

Our measurements were performed with the OPO cavity closely resonant to all three fields. 
We estimate that the OPO cavity detunings were smaller than 0.05 relative to its bandwidth 
for twin beams. In this situation, amplitude and phase quadratures decouple [19]. All analysis 
cavities were synchronously scanned, such that the same quadrature (amplitude, phase, or a 
linear combination of them) was acquired simultaneously for the three beams at each time. 

In our experimental investigation towards tripartite entanglement, we obtain data as presented 
in Fig. ID It shows the noise of the sum (full circles + line) and the difference (full line) of each 
pair of beams' quadrature fluctuations as functions of the analysis cavities' detunings A relative 
to their bandwidths (set to be equal during the experiment). The calibrated shot noise level 
is given by the dashed line. Fig. |2la) shows signal-idler combined noises. Fig. |2lb) considers 
signal-pump beams, and Fig. Wic) depicts idler-pump noises. All curves show the expected 
noise ellipse rotation: ampUtude quadrature noises are measured for large detunings (|A| > 3) 



and exactly on resonance (A = 0), while phase quadrature noises appear for detunings half- 
bandwidth off-resonance (A sa ±0.5). 

A third curve in each figure (open circles + line) shows the effect of including information 
about the third beam. They are obtained by subtracting the correction term from the sum curve, 
in the case of Fig. |2 a), or from the difference curve, in the cases of Figs.[2jb)-(c). The relevant 
regions in the corrected curves are located at A « ±0.5, corresponding to phase measurements 
[(a) a2^'+, (b) A^^^j and (c) A^^fjj]. All the curves presented in Fig. |2]come from a single 
measurement where all three fields' noises were recorded at the same time. 
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Fig. 2. Noise measurements of the terms appearing in the tripartite entanglement criteria 
of Ineqs. ([TJ-([3]l (pump power relative to threshold CT = P/Pth = 1.14). As functions of 
cavity detuning A, black full circles + line curves correspond to the noises in the sum of 
quadratures of two beams, solid red line to their difference, and blue open circles + line to 
a corrected noise, in which information from the remaining beam is included. Amplitude 
quadratures are measured at large cavity detuning |A| > 3 and at A = 0; phase quadratures 
are measured at A = ±0.5. (a) The considered pair of beams is signal and idler, and the 
corrected sum curve includes the term jSo coming from correlations with the pump beam; 
(b) signal-pump beams and dijference correction using idler; (c) idler-pump beams and 
dijference correction using signal. The corrected sum in (a) presents the first direct obser- 
vation of a quantum correlation among the phases of pump, signal, and idler: the noise is 
26% below the SQL. 

We begin our analysis by Fig. |2ja), devoted to signal-idler correlations and to the noise 
correction by pump beam. The well-known squeezing in /3_ is seen for large detunings in the 
difference curve, with noise value A^p_ = 0.45(1). The sum curve at A = ±0.5 gives the noise 



in equal to A^^+ = 1 .28 (2) . Although no squeezing is observed in q+ in this situation (pump 
power relative to threshold a = P/Pth = 1-14), entanglement is present, since A^/3_ + A^^+ = 
1.73(3) < 2, violating the Duan et al. inequality [16]. However, squeezing can be obtained 
in the three-field composite observable q'^, with the value ^q\ = 0.84(2), according to the 
corrected sum curve at A = ±0.5. This is the first direct observation of a quantum correlation 
among the phases of pump, signal, and idler Ineq. ([T]i then provides the expected violation of 
the first tripartite entanglement inequality, Vq — /^P- +A^^'^ = 1.29(2) < 2. 

Figs-Eb) and (c) present pump-twin correlations, and their difference noise corrected by the 
remaining twin. These two are very similar, as expected. The first term appearing in Ineq. (|2]i is 
measured in the sum curve of Fig.|2b). It is shot noise limited, A^/3oi = 1.03(3), as expected 
by the theoretical model, for this pump power The difference of their phases A^^oi presents 
excess noise, which is reduced by the correction term j32. The corrected difference phase noise 
^Qj, however, is not squeezed as expected, but shot noise limited as well, A^^qj ~ 1 .01 (4). This 
saturates Ineq. (|2]i, V\ = A^poi +^^^oi = 2.04(5) w 2. A similar result holds for Fig. |2tc). In 
this case, pQ2 presents a small amount of excess noise, A^po2 = 1.12(2), and is shot noise 
limited, A^^;,! 0.97(3), resulting in V2 = ^^PQ2+ ^^q'^i = 2.09(4) > 2 for Ineq. Q. 

4.1. Unexpected noise features 

Contrary to our expectations [13], we could not observe squeezing in A^^q^ or in A^^qj. We 
performed several measurements similar to those presented in Fig. |2] and recorded the values 
of noises for various pump powers. Each set of curves taken in one measurement scan, as the 
ones presented in Fig. |2] gives seven quantities of interest: A^/3_, A^^+, A^po> A^^O; 
A^;3oi, and A'^qj. Part of the results are presented in Fig. [3] where quadrature noises or their 
combinations are plotted as functions of pump power relative to threshold, a. 

In Fig. [3ja), the behaviors of terms belonging to Ineq. ([T]) are presented. The only noise 
curve that behaves in agreement with the theoretical model is A^/3_ (red circles). It's constant 
value is quantitatively compatible with what we expect from the other independently measured 
OPO parameters (mirror transmissions, total intracavity loss, OPO bandwidth, and detection 
efficiency [19]), A^/5_ ~ 0.46(1). Model prediction is the red solid line. 



Fig. 3. Noise behaviors as functions of pump power a. (a) Terms appearing in Ineq (TJ: 
the well-known squeezing in the subtraction of twins' intensity fluctuations A^j3_ is shown 
in the red circles, noise of twins' sum of phase fluctuations A^4+ are the black squares 
and the phase noise corrected by information from the pump A~^'^ are the blue triangles, 
(b) Amplitude (red circles) and phase (black squares) noise of the pump beam reflected 
by the OPO cavity. Solid lines are tentative fits to the ad hoc model explained in the text. 
Predictions differ from the usual ones only for the phase quadratures. Theoretical curves 
have the same colors as the experimental curve to which they relate. 
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Noise in the sum of twins' phases l^q+ (black squares) shows squeezing only very close to 
the oscillation threshold (CT < 1.1) and excess noise otherwise. The usual theoretical model pre- 
dicts that it should be squeezed for all values of a tested, reaching shot noise only for a ^ 1, if 
the incident pump beam is shot noise limited. The pump beam at the output of the filter cavity 
was verified to fulfill this requirement [11]. Yet, as we reported previously, excess noise is in- 
troduced in the intracavity pump beam by interaction within the non-linear crystal [11,14]. This 
was observed in the pump beam reflected from the OPO cavity, even when no parametric oscil- 
lation took place. The data presented in those papers was successfully modeled by introducing 
extra noise in the input pump beam. As we speculated in [11], this could be due to a non-linear 
refractive index effect, which would generate excess phase noise. The quasi-resonant cavity 
would then partially convert it into amplitude noise. In order to circumvent this effect, we de- 
cided to lower the pump finesse, which would minimize accumulated spurious phase shifts (the 
authors of ref. [9] claim this to be the main difference between their system, in which they do 
not observe this extra phase noise, and others). The OPO cavity finesse for the pump beam here 
is much lower than in our previous experiments. This modification entailed a lower reflected 
pump excess phase noise and completely suppressed the reflected pump amplitude noise but, 
alas, the mystery is far from solved, as described below. 

In Fig.|3lb), we present measurements of the reflected pump amplitude and phase quadratures 
during oscillation, respectively shown as red circles and black squares, as functions of C7. We 
could not account simultaneously for the data of Figs. [3] (a) and (b). We tried to mimic the 
spurious excess noise effect by introducing incident pump phase noise A^^j," in the equations, 
as an ad hoc free parameter Two different values of this parameter are needed in order to obtain 
the (already poor) agreement observed in Figs.[3](a) and (b). In order to obtain the black solid 
curve in Fig. |3] (a), we used the value A^q^, = 14; in order to do the same in Fig. [3] (b), we 
used A^q'p = 6. The reflected pump phase noise is lower than would be needed to account for 
the twins' excess phase noise. Furthermore, the reflected pump amplitude noise is even lower 
than expected by the linear model with no extra noise (already taking into account the finite 
detection efficiency). 

Next, we observe a correlation between the pump and twins' phase quadratures. In spite of 
its novel manifest quantum nature, we do not find agreement with the model. The blue solid 
line in Fig. [3] (a), was obtained by using ^^q'^j = 14. Since we are investigating a correlation 
between pump and twins' quadratures, and we could not find a single parameter to account for 
pump and twins separately, the poor agreement is not surprising. 

The noise features of the reflected pump phase quadrature are indeed more complicated than 
could be modeled by a single parameter. We measured this phase noise as a function of the anal- 
ysis frequency, for pump powers close to threshold and two different incident polarizations. In 
this measurement, the analysis cavity for the pump beam is locked to a detuning A ~ 0.5 and 
the photocurrents are analyzed by means of a spectrum analyzer (Resolution Bandwidth, RBW 
= 10 kHz, Video Bandwidth, VBW = 100 Hz). The phase matching conditions for the para- 
metric process are only fulfilled for one of these polarizations. We observed a rich spectrum, 
presented in Fig. HI with evenly spaced peaks [spacing approximately 150(10) kHz] of up to 
15 dB of excess noise, on top of a broad ~ 4 dB excess noise plateau. The peak structure is more 
pronounced for the "wrong" polarization, for which phase matching is not fulfilled [Fig.|4](b)]. 
At present, we cannot explain these noise features. Further investigation is needed, in order to 
reveal the physical mechanism behind such rich structure. 

4.2. New quantum correlations 

In spite of these unknown noise features, we could still check our experimental results against 
the tripartite criteria of Ineqs. (HJ - (O. We have already shown violation of Ineq. (HJ. Owing 
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Fig. 4. Reflected pump phase noise, as a function of thie analysis frequency. Pump power 
is close to threshold, but temperature is tuned to avoid triple resonance. In (a), the pump 
polarization is such that the phase matching conditions are fulfilled, whereas in (b) the 
orthogonal polarization is used. Measurements are normalized to the shot noise level (0 dB 
level). 
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Fig. 5. Behaviors of the noise terms appearing in Ineq (|2} as functions of pump power 
CT. The red circles correspond to the first term in the inequality and the blue triangles to 
the second one. The former presents the observation of amplitude quantum correlations 
between one of the twin beams and the reflected pump beam: there is ~ 12% of squeezing 
for higher a values. The solid lines are the predictions of the model explained in the text, 
without added noise. Theoretical curves have the same color as the experimental symbol to 
which they relate. 



to the similarity between the twin beams, it suffices to investigate Ineq. (|2]i. Our experimental 
results, as a function of CJ, are presented in Fig. [5] A new quantum correlation, among the 
amplitudes of the pump and the signal beams (manifest in A^/3oi, red circles), is observed for 
the first time in the above-threshold OPO. We obtain squeezing for higher values of pump 
power, A^poi ~ 0.88(2). Since it only involves amplitude quadratures, good agreement with 
the theoretical model (full red Une) is obtained (no excess noise is input in the model). The sum 
of this term with A^^qj (phase difference between pump and signal phases, corrected by idler 
phase noise) corresponds to the left hand side in Ineq. (|2|i. Data for A^^qj are presented as blue 
triangles in Fig. |5] The solid blue line corresponds to the theoretical model with no extra noise 



(adding extra noise does not provide good agreement). For a ~ 1 . 15, we approach the classical 
limit of Ineq. with the lowest value given by Vi =2. 04(10). 

Even by including spurious excess noise in the model, we would expect lower noise values 
for A^^Qj, leading to the predicted tripartite entanglement. Although we observed two new 
quantum correlations among the pump and the twin beams, genuine tripartite entanglement 
still remains to be demonstrated. Our results suggest that minor modifications to the setup 
may suffice. Technical improvements, such as a reduction of intracavity spurious losses and of 
detection losses, entail better squeezing and better detection efficiency. 

The ability to measure any given quadrature of each field and their combinations, together 
with the broad range of pump powers relative to threshold that we can access, allow for com- 
plete characterizations of the OPO's noise properties. In the near future, we expect to determine 
the excess noise source, to circumvent it, and reliably generate tripartite entanglement. 

5. Conclusion 

We presented a new and thorough investigation of quadrature noise of pump, signal, and idler 
fields in the above-threshold OPO. Although the OPO theory is well established, an important 
part of it had never been tested above the oscillation threshold. We revealed an unnoticed extra 
noise in the system and have further characterized several of its features. It seems safe to say 
that only phase noise is generated in the intracavity crystal. The amplitude noise observed in 
previous experiments was due to phase-to-amplitude conversion by the quasi-resonant cavity. 
This effective refractive index should be frequency dependent. Thus, by changing the pump 
frequency, it may be possible to minimize it. All this is still insufficient to pin down the ultimate 
physical origin of this noise. More research is needed to solve this intriguing puzzle. 

Notwithstanding the unwelcome extra noise, we could demonstrate new quantum corre- 
lations among pump and twin beams. The predicted quantum correlation among the phase 
quadratures, which is a direct consequence of energy conservation, was observed for the first 
time. When operating further above threshold, pump depletion becomes more evident as we 
observe the quantum correlation between the pump and each twin's amplitudes. By combining 
these results, we approached the limit for demonstrating tripartite entanglement: Vb=1.29(4), 
yi=2. 04(10), and V2=2.09(7), where a violation corresponds to Vj < 2. Hopefully, by improv- 
ing the squeezing and the detection efficiencies, we will beat this limit. 

The approach to directly generate multipartite entanglement, when compared to its genera- 
tion by means of individual squeezers and linear optics, has the potential advantage of entan- 
gling fields of different colors. This advantage is exploited in our setup and we expect it to be 
the starting point for multicolor quantum networks. 
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